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We study strongly correlated electrons on a kagome lattice at 1/6 (and 5/6) filling. They are described by an 
extended Hubbard Hamiltonian. We are concerned with the limit \t\ <^ V U with hopping amplitude t, 
nearest-neighbor repulsion V and on-site repulsion U. We derive an effective Hamiltonian and show, with the 
help of the Perron-Frobenius theorem, that the system is ferromagnetic at low temperatures. The robustness of 
ferromagnetism is discussed and extensions to other lattices are indicated. 
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Introduction. Ferromagnetism in solids or molecules can 
be of different origin. The most common is spin exchange be- 
tween electrons belonging to neighboring sites. Polarization 
of the spins and formation of a symmetric spin state reduces 
the effects of mutual Coulomb repulsions of the electrons due 
to the Pauli exclusion principle. The physics is the same as 
for intra-atomic Hund's rule coupling, which also plays a sig- 
nificant role in the theory of ferromagnetism. The key in- 
gredient of this mechanism is the potential energy of repul- 
sive electron-electron interactions minimized by a symmet- 
ric spin state which is better at keeping electrons apart. This 
should be contrasted with the standard superexchange mecha- 
nism for antiferromagnetism where the kinetic energy of elec- 
trons is optimized instead. Hence it is often the competition 
between the potential and kinetic energies that determines the 
"winner". This physics is illustrated, in its extreme limit, in 
the case of flat-band ferromagnetism 111 Mielke pointed 
out that electrons in a half-filled flat band become fully spin- 
polarized for any strength of the on-site repulsion U. (One 
could also think of this effect as an extreme case of the Stoner 
instability in metals.) 

Therefore it might appear surprising that ferromagnetism 
can also originate from purely kinetic effects. A prominent 
example is the ferromagnetic ground state (GS) discussed by 
Nagaoka jsl which is due to the motion of a single hole in 
an otherwise half-filled Hubbard system. The argument based 
on the application of the Perron-Frobenius theorem shall be 
presented later Although it is only valid in the limit of the 
infinite on-site Hubbard repulsion (to exclude the possibility 
of double occupancy) on a finite lattice, it demonstrates how 
ferromagnetism can result from the motion of the electrons or 
holes. The same theorem is also the basis of ferromagnetism 
due to three-particle ring exchange, a process first pointed out 
by Thouless [4] in the context of ^He (following the original 
observation by Herring jstl) and later also studied in the con- 
text of Wigner glass 0] and frustrated magnets |f^. In both 
cases, the ferromagnetic GS has the smoothest wavefunction 
and hence lowest kinetic energy. 

Our introduction would not be complete without mention- 
ing some other sources of feiTomagnetism such as the RKKY 
interaction in metals or double-exchange (e.g., in mangan- 
ites) to name a few. In this paper, however, we will be con- 



cerned with ferromagnetism of kinetic origin. In particular, 
we demonstrate that fermions on a partially filled kagome lat- 
tice which are described by an extended one-band Hubbard 
model in the strong correlation limit have a ferromagnetic GS. 
(The difference with Mielke's flat-band ferromagnetism is dis- 
cussed later in the paper) Again, the physics discussed here 
is motivated by the Perron-Frobenius theorem, but otherwise 
is quite different from Nagaoka's and Thouless' examples. 

Model Hamiltonian. We start from an extended one-band 
Hubbard model on a kagome lattice with on-site repulsion U 
and nearest-neighbor repulsion V. Using second quantized 
notation, the Hamiltonian is written as 

(1) 

Here the operators c^^ (cjg.) annihilate (create) an electron 
with spin a on site i. The density operators are given by 
Tii = nil + ^^'■h ^icr = cl^c-a- The notation refers 
to pairs of nearest neighbors. 

We first focus on the case of 1/6 filling (i.e., one electron 
per three sites). In the limit of strong correlations, when 
\t\ <^ V U and U ^ oo, the possibility of doubly oc- 
cupied sites is eliminated. First we assume that t = 0. In 
that case the GS is macroscopically degenerate. All config- 
urations with precisely one electron of arbitrary spin orienta- 
tion on each triangle are GSs (see Figs. [TJa) and (b)). It is 
helpful to consider the honeycomb lattice which connects the 
centers of triangles of the kagome lattice. Different GS con- 
figurations on the kagome lattice correspond to different two- 
colored (spin) dimer configurations on the honeycomb lattice 
(particles are sitting here on links, see Fig.[T](c) and (d)). They 
are orthogonal because any wavefunction overlap is neglected. 

When t 0, this GS degeneracy is lifted. In the lowest 
non-vanishing order in t/V, the effective Hamiltonian acting 
within the low-energy manifold spanned by the states with no 
double occupancy and exactly one electron per triangle be- 
comes 

^hex = -5 E (lQ)(C>l + lO)(Dl+H.c.) 

(2) 
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FIG. 1: Panels (a) and (b) show two different configurations satisfy- 
ing the constraint of zero or one electron per site and one electron of 
arbitrary spin per triangle. The arrows indicate possible ring-hopping 
processes. An equivalent colored dimer representation on a honey- 
comb lattice is shown in (c) and (d); the colors encode particle spins. 
Ring exchanges conserve the parity of the number of dimers on the 
sublattice shown in panel (e). Panel (f) shows the 24-site cluster used 
for the exact diagonalization. The dimers are arranged to maximize 
the next-nearest neighbor spin interaction along the dashed lines. 



with g = Qt^ /V^ . Here the Hamiltonian is written in terms of 
dimers on a honeycomb lattice and the sum is performed over 
all hexagons. The sum over the three symbols is taken over 
all possible color (spin) combinations of a flippable hexagon. 
Particles hop either clockwise or counter-clockwise around 
the hexagons. These processes can lead to different config- 
urations, depending on the colors (spins) of the dimers. We 
observe that i/hex does not cause a fermionic sign problem. In 
particular, the local constraint of having one fermionic dimer 
attached to each site allows for an enumeration of dimers such 
that only an even number of fermionic operators has to be ex- 
changed when the matrix elements of iJhex are calculated. 

If one were to ignore the spin degrees of freedom (the col- 
ors of the dimers), the model would be equivalent to the quan- 
tum dimer model (QDM) studied in Ref. fl] . Similarly to the 
QDM on a square lattice |9], the effective Hamiltonian ^ 
conserves certain quantities - winding numbers - and con- 
nects configurations only when they belong to the same topo- 
logical sector (For the case of periodic boundary conditions, 
the winding numbers are defined by first orienting dimers so 
that the arrows point from the A to B sublattice, and second, 
by counting the net flow of these arrows across two indepen- 
dent essential cycles formed by the dual bonds.) The GS of the 



QDM was found to be three-fold degenerate in the thermody- 
namic limit, corresponding to the valence bond solid (VBS) 
plaquette phase with broken translational invariance. In what 
follows, we shall investigate the effects of quantum dynamics 
- the ring-exchange hopping of electrons (dimers) - on spin 
correlations. Note that i/hex has no explicit spin dependency 
and conserves both S^^^ and the total spin Stot- 

Ferroniagnetism from the Perron-Frobenius Theorem. In 
short, the Perron-Frobenius theorem states that the largest 
eigenvalue of a symmetric nxn matrix with only positive el- 
ements is positive and non-degenerate, while the correspond- 
ing eigenvector is "nodeless", i.e., can be chosen to have only 
positive components. (For a simple proof of this theorem, see 
e.g. ifioll .) Applying this theorem to the (finite-dimensional) 
matrix exp {—tH) (for any r > 0), one concludes that if all 
off-diagonal matrix elements of the Hamiltonian H are non- 
positive and the Hilbert space is connected by the quantum 
dynamics (meaning that any state can be reached from any 
other state by a repeated application of H), then the GS is 
unique and nodeless. It is important to remember that the the- 
orem only works for systems with a finite-dimensional Hilbert 
space. 

To show the relation of this theorem to ferromagnetism, we 
now sketch the argument for Nagaoka's ferromagnetism in the 
GS of an infinite U Hubbard model (Eq. © with V = Q, 
U oo) with a single mobile hole (after Refs. Ill ll 11211 ). 



Denote a state with a single electron or hole on a site i 
as |i, a) where a = {ai, . . . ai-i,ai+i, . . . aN} is the spin 
configuration of electrons. We use the convention \i, a) = 
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^^|0). (No double 



occupancy is allowed if [/ — > oo.) In this basis, the matrix el- 
ements of the hopping Hamiltonian are either t, for the states 
related by a single hop of the hole between two neighbor sites, 
or otherwise. The Hamiltonian commutes with both S^^^ and 
S'fgj. Our chosen basis consists of eigenstates of S^^^ but not 
Stot' hence we can immediately separate the Hilbert space into 
sectors of fixed S^^^. Within each sector, the Hamiltonian ma- 
trix has exactly z (the coordination number) nonzero entries in 
each row and each column. A direct inspection shows that a 
vector whose entries are all 1 is an eigenvector with the eigen- 
value zt. If t < and tunneling of a single hole satisfies the 
connectivity condition, the Perron-Frobenius theorem applies 
and hence such a state is the GS (there can be no other state 
with positive coefficients only that is orthogonal to this one), 
which is unique for a finite system. (We remark that the sign 
of t can always be changed on a bipartite lattice.) Clearly, this 
state is fully spin-polarized in the S'j^, = S'^ax = ±iV/2 sec- 
tors, and since the Hamiltonian commutes with S^^^, the state 
with = (S'niax+l)<S'max must havc the same energy in every 
5tot sector But we already saw that the states with the energy 
£ ~ zt are unique GSs in every sector, hence they must have 
the same S^^^ = (^max + l)S'max, QED. The obvious pitfalls 
may come from taking the thermodynamic limit or violating 
the connectivity condition: in both cases the nodeless, fully 
spin polarized state remain a GS but no claims can be made 
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about other potential GSs. 

Turning to our case, we remark that the sign of the plaque- 
tte flip in Eq. ^ can be always chosen negative, irrespective 
of the sign of the original tunneling amplitude t - this is just 
a matter of a simple local gauge transformation j^]. Specifi- 
cally, the sign of g in Eq. Q can be changed by multiplying 
all configurations C with the color-independent factor i'^^'^^ 
where !^(C) is the number of dimers on the sublattice shown 
in Fig. [T](e). This fact might appear surprising though since 
the actual sign of t can typically be gauged away only for 
the cases of bipartite or half-filled lattices. The reason the 
sign of t turns out to be inconsequential in our case of the 
(non-bipartite) kagome lattice away from half-filling is due 
to the constrained nature of the ring exchange quantum dy- 
namics of Eq. (|2]). We therefore choose all off-diagonal ma- 
trix elements of i^hex to be non-positive. This by itself is not 
yet sufficient to apply the Perron-Frobenius theorem since the 
quantum dynamics of dimers on a (bipartite) honeycomb lat- 
tice is explicitly non-ergodic: as we have mentioned earlier, 
the Hilbert space is broken into sectors corresponding to the 
winding numbers which are conserved under any local ring 
exchanges. On the other hand, the ring-exchange dynamics 
of dimers given by Eq. (|2]i is ergodic within each sector ifisll . 
Therefore we consider each topological sector separately. The 
argument is very similar to the one presented earlier for Na- 
gaoka's ferromagnetism. For the S^^^ ~ Smax — ±-/Ve/2 spin 
sectors, the GS is unique, fully spin-polarized, and all ele- 
ments of its eigenvector are positive. 

For all other S^^^ sectors, however, the situation appears 
more complicated at first sight. The reason is a much big- 
ger configuration space - essentially, we are now dealing with 
two-color dimer configurations. A given state can now be 
connected by a ring-exchange Hamiltonian to a larger num- 
ber of states than it would if all dimers had the same color 
(spin). We formalize this by introducing the notion of de- 
scendant states |Cf ), k = 1 . . .2^' - two-color dimer con- 
figurations obtained from the uncolored "parent" configura- 
tion \Ci) by simply coloring its dimers (i.e., assigning spins). 
The subspace of descendant states can be partitioned accord- 
ing the conserved S^^^. The resulting sectors, in general, have 
different dimensionality D{S^^^) equal to the number of dis- 
tinct permutations of spins (colors). A crucial observation 
is that Efe(^^f l^hexICf ) = (C,|i/hex|C,) for any ^, j, m. 
The immediate consequence is that if j^'g) = |^o(-^e/2)) = 
J2i 7i|C'i} is the GS in the S^^^ = S'^ax = ±-^e/2 spin sector, 
then l^oiS^o,)) = ^-^/'(^fot)E.7.ELl^^^) is an eigen- 
state with the same energy in any other S^^^ spin sector (The 
sum over descendants k is performed only within a given spin 
sector) The SU(2) symmetry of the effective Hamiltonian ^ 
once again implies that l^'o('S'tot)) is a GS and is fully spin 
polarized. 

Notice that the uniqueness of such a GS relies on the ergod- 
icity of the Hamiltonian within each S^^^ spin sector Numer- 
ical studies on finite clusters up to 48 kagome lattice sites in- 
cluding different geometries show that this is in fact the case. 
Unfortunately, we were not able to provide a rigorous analyt- 



ical argument. Should it turn not to be the case, it would open 
a possibility for degenerate GSs in S^^^ ^ S'max spin sectors. 
Still, at least one of the GSs is fully spin polarized fl^ . 

While the Perron-Frobenius argument applies to finite sys- 
tems, it does not withstand the thermodynamic limit. In par- 
ticular, it is known that in the thermodynamic limit the ground 
state of the QDM is in the three-fold degenerate plaquette 
phase |@]. We suspect that ferromagnetism survives this limit 
and coexists with such a broken symmetry state; a conclusive 
resolution of this point remains a subject of further research. 

The ferromagnetic GS which we find here should not be 
mistaken with Mielke's flat -band ferromagnetism |[ll|2l- In 
fact, Mielke has shown that a positive-?/ Hubbard model with 
V = on a kagome lattice at 5/6 filling has a fully spin- 
polarized GS. A detailed discussion of the differences to our 
case can be found below. 

Stability of kinetic ferromagnetism. In order to test the ro- 
bustness of the ferromagnetic GS, we introduce by hand an 
additional next-nearest neighbor interaction, in the spirit of 

H' = i?hex + J i^'^i - ■ 

By adding such a term, we attempt to frustrate the ferromag- 
netic state. Indeed, this term favors configurations such as 
the one shown in Fig. [Tff) which maximize the spin interac- 
tions by having electrons on the same sublattice. Not only 
the resulting charge order is expected to suppress the kinetic 
mechanism for ferromagnetism, the ferromagnetic order itself 
is now suppressed by antiferromagnetic fluctuations favored 
by the spin interactions for J > O.By gradually increasing 
J /g towards strong antiferromagnetic coupling, we can esti- 
mate the stability of the ferromagnetic GS under the presence 
of short ranged perturbations. 

Numerical Results. We calculate by means of exact diago- 
nalization the GS of a two-color dimer model on a 18- and 24- 
site honeycomb lattice. The system corresponds to a 1/6 filled 
kagome cluster with 27 and 36 sites, respectively. The cal- 
culations on the clusters of different size and geometry show 
qualitatively the same results. 

The ground-state energies of the different Sz sectors are 
degenerate as long as J / g < J/g)c ~ 0.2, as shown in 
Fig. 12] This demonstrates the robustness of ferromagnetism 
induced by ring-hopping processes. Above the transition 
point [J / g)c ~ 0.2 antiferromagnetic spin fluctuations are no 
longer suppressed. The ground-state degeneracy of the differ- 
ent Sz sectors is lifted and the true GS is the one with the low- 
est \Sz\- The gain in kinetic energy decreases correspondingly 
since the spin fluctuations cause nodes in the wavefunction. 
The expectation value of i/spin shows a jump at {J/g)c (see 
Fig. |2|. Note that a small second jump at J/g « 0.45 found 
for the 18 site cluster is an effect of geometry. The observed 
findings demonstrate a considerable robustness of ferromag- 
netism generated by kinetic processes. In the limit J/g oo, 
the kinetic processes are unimportant and the GS is that of a 
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FIG. 2: Exact diagonalization of the two-color dimer model on a 18- 
site (left panels) and 24-site (right panels) honeycomb cluster. The 
upper panels show the GS energies of different sectors as a func- 
tion of next-nearest neighbor coupling J/g. The lower ones show 
the GS expectation values of the spin part of the Hamiltonian. 



Heisenberg antiferromagnet on a kagome lattice. One of the 
configurations which maximize the spin interactions is shown 
in Fig. [11(f). 

The above considerations can also be applied to the case of 
5/6 filling due to the arbitrary choice of sign of g in Eq. Q. In 
addition, a ferromagnetic GS is found numerically for a fill- 
ing factor of 1/3. With two occupied sites per triangle, we 
obtain a fully packed loop covering of the honeycomb lat- 
tice instead of a dimer covering. As before, there is no sign 
problem in the strong coupling limit. However, the effective 
Hamiltonian is not ergodic in the 1/3 filled case and thus the 
Perron-Frobenius theorem does not rule out other GSs which 
are not fully spin polarized. Numerical studies confirm a fer- 
romagnetic GS which is much less robust - antiferromagnetic 
fluctuations occur already at very small ratios {J/g). A more 
detailed discussion is left to an extended version of this paper 

We conclude by reiterating the difference between the two 
mechanisms for ferromagnetism in a Hubbard model on the 
kagome lattice: the one presented here and the flat-band 
mechanism discussed in Refs. The flat-band fer- 

romagnetism was demonstrated for the case of y = in the 
Hamiltonian ([l]), while our mechanism requires V ^ oo. 
In Mielke's case, ferromagnetism has been predicted for the 
range of fillings between 5/6 and 11/12 and any value of 
U > {). The connection between the sign of the tunneling 
amplitude t and the electron concentration is crucial. This 
is because the tight binding model on a kagome lattice has 
one of its three bands completely flat: the lowest band for the 



case of t < or the highest band for the case of t > 0. On 
the other hand, the mechanism presented here is insensitive to 
the sign of t and, as we have already mentioned in the intro- 
duction, is of kinetic rather than of potential origin. Further- 
more, Perron-Frobenius theorem is not applicable in Mielke's 
case lll4ll . instead the proof was based on graph-theoretical 
methods. By contrast, the kinetic ferromagnetism studied in 
this letter relies crucially on strong electron-electron repul- 
sion t/ ^ y ^ cxd; hence it belongs to a different class from 
flat-band ferromagnetism. Despite certain similarities, it must 
also be distinguished from Thouless' three-particle ring ex- 
change mechanism: the crucial difference is that a standard 
three-particle ring exchange leaves the particles at the original 
locations, it simply cyclically permutes them. In our case the 
particles actually move; the initial and final configurations are 
distinct. 

Thus the mechanism found in this letter represents a new 
generic type of kinetic ferromagnetism. An interesting ques- 
tion is to what extent the ferromagnetism of the form found 
here can be extended to other lattice structures. A particularly 
interesting case is the pyrochlore lattice. For the strong cor- 
relation limit at 1/8 filling one can use the PeiTon-Frobenius- 
based argument as well, implying a ferromagnetic GS. Again, 
this and other cases will be discussed in the extended version 
of this paper. 
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